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Summary. I present some new results, concerning identification of con- 
volution kernels in integrodifferential equations of parabolic type. I consider, 
in particular, fully nonlinear equations with kernels depending only on time 
and linear equations with kernels which can depend also on some of the 
space variables. 


Riassunto. Presento alcuni risultati nuovi sull’identificazione di nuclei 
di convoluzione in equazioni integrodifferenziali di tipo parabolico. Con- 
sidero, in particolare, equazioni totalmente non lineari con nuclei che dipen- 
dono solo dalla variabile temporale ed equazioni lineari con nuclei che pos- 
sono dipendere anche da qualcuna delle variabili spaziali. 
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Introduzione I problemi inversi consistono nello studio di equazioni 
di vario tipo (tipicamente differenziali o integrodifferenziali), in cui ciò che 
tradizionalmente viene assunto come dato é conosciuto solo parzialmente. 
Questa carenza di informazione é compensata dalla conoscenza di qualche 
aspetto della soluzione, tradizionalmente assunta come totalmente incognita. 

In [9] considerai vari esempi di problemi differenziali e integrodifferen- 
ziali. Qui mi occuperé esclusivamente di equazioni integrodifferenziali di 
tipo parabolico ed esporrò qualche nuovo risultato, ottenuto in collabo- 
razione con F. Colombo, del Politecnico di Milano, e con A. Lorenzi, dell’ 
Università di Milano. Più specificamente, mi occuperò della determinazione 
(contemporaneamente alla soluzione) di nuclei di integrazione per equazioni 
integrodifferenziali di tipo parabolico. L’esempio pit semplice del tipo di 
problema che vado a considerare é il seguente: consideriamo l’andamento 
della temperatura u in un corpo {2 costituito da un materiale omogeneo con 
memoria. Se u é mantenuta costante sul bordo 09 di Q (per semplicità, 
supponiamo u(t,r) = 0 per ogni t € Re per ogni x € 0N) e se é noto il 
suo andamento prima di t = 0, un modello matematico standard di u (con 
qualche semplificazione) é costituito dal problema 


du(t,0) = Au(t,2)+/h(t- s)Au(s,r)ds + f(t,r) (t,7) [0,T]xQ, 
u(0, x) = uo(7), re, 

u(t,z)=0, t € [0,T],r € 09, 

(1) 

ove h :]0,7] — R rende conto del fatto che il materiale é con memoria. 

Assegnata hh, l’equazione (1) nell’incognita u é stata abbondantemente stu- 

diata. Nella pratica, però, h é raramente conosciuta. A compensare questa 

carenza di informazione, é spesso possibile determinare l’ andamento di un 

certo funzionale @ in corrispondenza della soluzione u. Per esempio, é tal- 

volta possibile fissare un certo punto 3 € {2 e osservare l'andamento della 

temperatura u(t, r1) al variare di t in [0,7]. In tale caso, l’indeterminatezza 

di — é compensata dalla conoscenza di u(t, 71) per ogni t € [0,7]. Più in 
generale, supponiamo assegnato 


H(u(t, .)) = g(2), (2) 


ove è é un opportuno funzionale. Una versione astratta del problema (1) fu 
studiata nel lavoro [10]. 

Un problema totalmente non lineare Consideriamo ora una gene- 
ralizzazione astratta totalmente non lineare del problema (1). Il problema é 
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il seguente: determinare 7 > 0 e due funzioni u : [0,7] + D, h: [0,7] + R, 
tali che 


d,u(t) = G(u(t)) + (o A(t — s)F(u(t),u(s))ds + f(t), t € [0,7], 
u(0) = vo, (3) 
P(u(t)) = g(t), te [0,7]. 


Qui D e X sono due spazi di Banach, D C X con immersione contunua, 
G:D-X,F:DxD- X,f:[0,7] > X, un € Dsi suppongono 
noti. L'ignoranza del nucleo h é compensata dalla conoscenza di d(u(t))(= 
g(t)), con ® funzionale lineare e continuo su X. Il problema (2) si suppone 
parabolico nel seguente senso: che il differenziale secondo Frechét G'(uo) di 
G in corrispondenza di up (che é un operatore lineare e continuo da D a X) 
é un operatore settoriale in X. Ciò significa quanto segue: 

(I) l'insieme risolvente p(G'(uo)) di G'(uo) contiene {A € C : |Aj > 
R, |ArgA| < 5}, per qualche R > 0; 

(II) esiste M > 0 tale che, per ogni À nella regione suddetta, ||[(A — 
G'(uo)) lex) < MAI. 

In un lavoro sostanzialmente ultimato in collaborazione con F. Colombo 
(vedi [1]) abbiamo studiato il problema (2) utilizzando metodi di regolarità 
massimale in spazi di Sobolev di ordine non intero. Tali spazi sono preferi- 
bili ai più classici spazi di Sobolev di ordine intero in quanto ammettono 
una teoria della regolarità massimale assai più semplice (vedi, ad esempio, 
[6] per risultati di regolarità massimale in spazi di Sobolev di ordine intero). 
Precisiamo che per risultato di regolarità massimale intendiamo un risul- 
tato di isomorfismo vettoriale e topologico. Per enunciare il risultato che ci 
interessa, é necessario introdurre alcune definizioni preliminari. 

Cominciamo col richiamare la definizione di spazio di Sobolev di ordine 
B €]0, i[. Sia p € [1, +00], sia X uno spazio di Banach e T > 0. Poniamo 


IP(O.T-X): (T( pt MOI a 
8 {f e LP(0, 5A): So (So Tiegrrant ds < +00) 
W°P(0,T;X):= se p< +00, (4) 
: C8(0,T; X) se p= +00. 


Qui Cf(0, T; X) indica il classico spazio delle funzioni hòlderiane di ordine 8. 
Sia A un operatore settoriale in X. È noto che A é il generatore infinitesimale 
di un semigruppo (7(t)):>o non necessariamente fortemente continuo in 0 
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(non supponiamo che il dominio di A sia denso in X). Sia # €]0, 1[, p € 
[1, +00]. Poniamo 


Da(B,p):= {re X:t+ |!" AT(t)e]|x € L?(0,1)}. (5) 


Possiamo ora enunciare il risultato di regolarità massimale che utilizze- 
remo per il problema inverso: 


Teorema 1 Sia A un operatore settoriale nello spazio di Banach X. Con- 
sideriamo il problema di Cauchy 


| u'(t) = Au(t)+f(t), te [0,T], 
u(0) = uo. 


(6) 


Le seguenti condizioni sono necessarie e sufficienti affinché il problema (6) 
abbia una soluzione (unica) u appartenente a W!+8P(0,T; X)NW?:P(0, T; D 
(A)), per assegnati p €]1, +00) e f €]0, 1\{}}: 

(I) f € WBP(0,T;X); 

(II) 

Da(B+1-3,p) se B<l, 
mula D(A) se p>Ì 

(III) se B> }, allora Au + f(0) € Da(6— 3,p). 

[1] contiene una dimostrazione diretta del teorema 1. Una diversa di- 
mostrazione si può anche ottenere da certi recenti risultati di G. Di Blasio 
(vedi [5]). 

Torniamo ora al problema (3). La difficoltà principale sta nel fatto che la 
prima condizione é un’equazione di Volterra di prima specie nell’incognita 
h. Ora, una strategia standard per passare da un’equazione di Volterra 
di prima specie a un’equazione di seconda specie é derivare l’equazione. 
Affinché ciò sia possibile, si é costretti a cercare una soluzione vw assai 
regolare. Pit precisamente, cercheremo una soluzione (u,h) di (3) con 
u € W2+4P(0,7;X)M W!+4P(0,7; D) e ha € WSP(0,7) (a valori scalari), 
per qualche # €]0, 1[, p €]1, +00], 7 > 0. Vale il seguente 


Teorema 2 Consideriamo il problema (3). Fissiamo preliminarmente p € 
]1, +00] e 0 €], 1\{}} Introduciamo le seguenti ipotesi: 

(a) G € C*(D; X), F € C3(D x D; X); 

(b) uo € D; 
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(c) A := G'(uo) é un operatore settoriale in X; 
(d) f e W!+42(0,T; X); 
(e) 


sil 
cao) +10) ef PAP*) pr) se f< 


(f) ge W?+AP(0,T); 

(9) ®(uo) = g(0), ®(G(uo) + f(0)) = 9/(0); 
(h) x := ®(F(uo,vo)) # 0; 

(i) se BP > 3 e 


H:= x 1{9"(0) -— B[A(G(uo) + f(0)) + f'(0)], 
allora 


A[G(uo) + £(0)] + f"(0) + HF (uo, uo) € Da(B— Lp) 


(1) F(uo, uo) € D. 

Allora il problema (3) possiede almeno una soluzione (u, h) € (W?+4P(0, 
T;X)NW!+4P(0,7; D)) x WAP(0, T), per qualche  < T. Inoltre, se F(u(t), 
uo) € D per ogni t € [0, 7], tale soluzione é anche unica in (0, 7]. 


Osservazione Si verifica facilmente che le ipotesi (b), (c), (f) e (9) 
sono necessarie per ottenere una soluzione con la regolarità voluta. Sotto la 
condizione (h), si verifica poi che anche le ipotesi (e) e (i) lo sono. Cié segue 
abbastanza facilmente dal teorema 1 (vedi la dimostrazione successiva). 


Cenno della dimostrazione del teorema 2 Ammettiamo che il pro- 
blema (3) possieda una soluzione (u, h) appartenente a (W?+4P(0,7;X)n 
W1+5P(0, 7; D))x WAP(0, 7), per qualche 7 < T. Poniamo v := du. Allora, 
derivando la prima equazione e ricordando che A := G'(uo), si ottiene che 
(v, h) € (W!+8(0,7;X)N WAP(0,7; D)) x W8-P(0, 7) e risolve il problema 


O v(t) = Av(t) + h(t)F(uo, uo) + f'(t) +P(v,h)(t) te [0,1], 
v(0) = G(uo) + f(0), (7) 
d(v(t)) = 9'(t), te [0,7], 


ove 


P(v, h)(t) := [G"(uo+1+v(t)) — AJu(t) + h(t)[F(u0+1+v(t), uo)+ 
—F(uo,uo)] + /{ A(t— s)dF(uo +1+U(t),uo+14 v(s))(v(t), v(s))ds. 
(8) 
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Qui abbiamo indicato con + la convoluzione su [0, +00[. Si osservi che u = 
uo +1+v. La prima equazione in (7) é stata scritta in modo da evidenziare il 
fatto che pensiamo al secondo membro come a una perturbazione di Av(t) + 
h(t)F(uo, uo) + f'(t). i 

Applichiamo ora ® alla prima equazione. Usando l’ultima condizione in 
(7) e l'ipotesi (kh), otteniamo il sistema 


d;v(t) = Av(t) + h(t)F(uo, uo) + f‘(t) +P(0,h)(t) t € [0,7], 
v(0) = G(uo) + f(0), (9) 
h(t)= x tg"(t) — D(Av(1) + P(v,h)(t) + f"(2))], t € [0,7], 


Si tratta ora di risolvere il sistema (9). La formula di variazione delle 
costanti consente di riscrivere la prima condizione in (9) nella forma 


v(t) = T(t)[G(uo) + f(0)] + Jo T(t — s)f"(s)ds+ 
+ SE T(t— s)[h(s)F (uo, vo) + P(v, h)(s)]ds. 


Ora, l’ipotesi (/) consente di considerare, almeno nel caso 8 < 1, anche 
h(t)F(uo, uo) come una perturbazione, in quanto vale il seguente 


(10) 


Lemma 1 Siano 6 < i, he WfP(0,7) eye D. Allora 


If TC- Ud wsrtorD) £ IA Ilwe107) 
con limm(7) = 0. 


Se i termini di perturbazione fossero nulli, si potrebbe utilizzare il teo- 
rema 1, determinare una soluzione v € W!+8P(0,7;X)N W£P(0, 7; D) del 
sistema costituito dalle prime due equazioni e ricavare A dall’ultima con- 
dizione in (9). Le ipotesi sui dati iniziali (d), (e) e (i) sono proprio ciò che 
serve per poter applicare il teorema 1. Pi in generale un metodo pertur- 
bativo standard, basato sul teorema delle contrazioni, permette in ogni caso 
di trovare una soluzione su un intervallo [0, 7] con 7 ” piccolo”. O 


Un’ applicazione del teorema 2 Vediamo ora un'applicazione del teo- 
rema 2 a una generalizzazione di un problema di dinamica delle popolazioni 
(vedi [7]). Consideriamo il seguente problema misto quasi lineare: 
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Au(t,x) = d(u(t,x))Azu(t,x)) + bu(t, 2) Sh(t- s)u(s, r)ds + f(u(t,2)) 
te [0,7], re, 
u(t,r)=0,t € [0,7], r € 90, 
u(0,r) = uo(7), 7 EQ, 
Ja d(2)u(t, 2)dr = g(t), t € [0, 1]. 

(11) 
nelle incognite u : [0,7]xN + Re h:{0,7] + R. Qui Né un aperto limitato 
e regolare in R" e db € R.\ {0}. Osserviamo ancora che l’indeterminatezza 
di h é compensata dalla prescrizione di SJa9(£)u(t, )dr. Supponiamo anche 
‘che: 

(h11) d € C*(R) e d(u)> 0 per ogni u in R; 
(h12) è € IP (9), conp>1V5; 
(h13) f € C*(R). 


Poniamo 
X:= D(Q), (12) 
D:=W?P(Q)NWP(0), (13) 
G:D- X,G(u) := d(u)Au+ f(u)vu € D, (14) 
F:DxD- X,F(uz,u2):= bu;ugVu;, ug € D. (15) 


Osserviamo che da (h12) segue che D é immerso con continuità in C(9) 
ed é facile verificare che G € C*(D; X), mentre F € C3(D x D; X). Se 
uo € D, si ha inoltre 


Av = G'(uo)v = d(uo)Av + [Auod'(uo) + f(uo)]v Vu € D. (16) 


È ben noto (vedi, ad esempio, [12]) che A, di dominio D, é un operatore 
settoriale in X. Si osservi che in questo caso si ha D= X. 

L'ultimo ingrediente che manca per poter applicare il teorema 2, é la 
specificazione degli spazi D4(0,p), con 9 €]0, 1[. Nel caso specifico, si può 
verificare che, con equivalenza di norme, per 0 £ dà 

20,9) = | B?0(9) se 0<0<5%; 


{u € B2.(9) ‘uan =0 se x «zl 19) 
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ove, con la scrittura Bi a(9), indichiamo lo spazio di Besov di indici cor- 
rispondenti. 

Per una dimostrazione di questo risultato, dovuto a P. Grisvard, si veda 
anche [8]. 

È a questo punto un facile esercizio, che lasciamo volentieri al lettore, 
applicare il teorema 2 allo studio del problema (11). 


Nuclei dipendenti anche da variabili spaziali e soluzioni globali 
Come qualcuno avrà osservato, il risultato astratto contenuto nel precedente 
teorema 2 é di tipo locale nella variabile temporale: in generale si é in 
grado di provare l’esistenza e unicità di una soluzione solo in brevi intervalli 
temporali. Si osservi che anche nel caso lineare la questione é tutt'altro che 
banale,in quanto, anche se si parte da un problema lineare in u, il problema 
inverso é non lineare nell’incognita (u, h) (si consideri, ad esempio, (1)). 

Nel caso lineare, con nucleo dipendente esclusivamente da t, l’esistenza di 
una soluzione globale é stata studiata in [11]. Un altro problema interessante 
é costituito dalla determinazione del nucleo ‘ nel caso in cui esso non dipenda 
esclusivamente dal tempo. 

Un risultato in questa direzione é stato ottenuto in [4]. Un altro risultato, 
di tipo astratto, é stato ottenuto in [3] e ha il pregio di garantire anche 
l’esistenza e l'unicità di una soluzione globale. Per darne un'idea, enuncio 
una semplice applicazione del risultato principale contenuto in [3]. 

Sia O il quadrato ]0, 1[x]0, 1 (C R?) e sia T > 0. Vogliamo determinare 
due funzioni U : [0,7] x T-+ Rek:[0,7] x [0,1] + R tali che 


aU(t,r,y) = Ault,r,y) + & k(t— s,x)Au(s,r,y)ds + f(t,2,y), 
te [0,7], (2,y) EQ, 
u(t,x,y)=0,t € [0,T], (c,y) € 99, (18) 
u(0, 2, y) = uo(£, y), (7,7) EI, 
ST P(c,y)u(t, 2, y)dy = g(t,2),t € (0, 7], e € (0,1). 


Qui abbiamo indicato con A l’operatore di Laplace nelle variabili (2, y). 
Osserviamo che stiamo cercando di determinare un nucleo k dipendente 
dalle variabili t e x. 

Indichiamo ora con A l’operatore definito come segue: 


D(A) := {ue H?(Q)N H}(A)NC(A) : Au e C(A)}, (19) 


Au:= Au Vue D(A). (20) 
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É noto che A é un operatore settoriale in C(Q) (vedi [2]) e D(A) C 
NesoCt (I), ove abbiamo posto 


CHA) := {fue C.D) : upon = 0). (21) 


Teorema 3 Sia 8 €]0, 1[\{}}. Supponiamo che siano soddisfatte le seguenti 
condizioni: 

(a) $ € CI); 

(5) f € CIH#([0,T);C@)); 

(6) uo € D(A), vo := Aug + Î(0)(= Auo + f(0)) € D(A); 

(c) So #(x, y)Auo(x,y)dy #0 Ve € [0, 1}; 

(d) g € C°+2([0, T]; C([0, 1})nC!+4([0, T]; C2([0, 1)), 9(t,0) = g(t,1)=0 
Vi € [0,7]; 

(e) Sa (2, y)uo(x,y)dy = 9(0,2), fl d(2, y)(Avo(2,9) + f(0,2,9))dy = 
dg(0, x) Ve € [0,1]; 

(#7) Avo +0.f(0) + koAuo € C?A(T), ove 


ko(x) == 06 9(0, x) Di ki d(2, y)A?uo(z, y)dy 
h d(7, y)Auo(x, y)dy 


Allora, il problema (18) possiede un'unica soluzione (u,k), con u € 
C?+2([0, T]; C(A) N CI+4([0, T]; D(A)) e k € CP([0, 7]; C((0,1))). 
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